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Travelling- wave solutions of the standard and compound form of Korteweg- 
de Vries-Burgers equations are found using factorizations of the correspond- 
ing reduced ordinary differential equations. The procedure leads to solutions 
of Bernoulli equations of nonlinearity 3/2 and 2 (Riccati), respectively. In- 
troducing the initial conditions through an imaginary phase in the travelling 
coordinate, we obtain all the solutions previously reported, some of them 

\& \ being corrected here, and showing, at the same time, the presence of inter- 

esting details of these solitary waves that have been overlooked before this 

J^l, . investigation. 
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1. INTRODUCTION 

X 

5_i ■ In this research we will study nonlinear Korteweg-de Vries-Burgers (KdVB) 

equations which play an important role both in physics and in applied math- 
ematics. First, we will consider the original KdVB equation 01212], which 
can be written as 

du d 3 u d 2 u du 
dt dx 3 dx 2 dx 
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Second, we will focus on the so-called compound KdVB equation that has 
an additional non-linear term [3J [SJ H3 [7] 

du d 3 u d 2 u du 2 du 

dt dx 3 dx 2 dx dx 

In both cases the parameters s, fx, a, (5 are real constants which take into 
account different effects, as non-linearity, viscosity, turbulence, dispersion 
or dissipation. As it is well known, these equations have been used as math- 
ematical models for the propagation of waves on elastic tubes [5], and in 
particular, the KdVB equation has been explicitly derived in the case of 
ion-acoustic shock waves in multi-electron temperature collisional plasma 
[2] . These ion-acoustic shock waves have been observed in an unmagnetized 
plasma ^Uj. Standard applications are encountered in hydrodynamics and 
in the theory of diffusive cosmic-ray shocks |11[ I12j . Other interesting and 
more recent applications refer to the dynamics of high-amplitude picosec- 
ond strain pulses in sapphire single crystals in which case a KdVB model 
has been shown to fit the observed behavior at all strains and temperatures 
under study ^3] an d explaining the temporal self-shift of small-amplitude 
dark solitons due to self- induced Raman scattering in optical fibers ^3] . 

We are interested in travelling- wave solutions of the Eqs. (fl | )— ([% ]) . that 
is, solutions of the form 

u(x,t) = 4>(0, i = x-vt, (3) 

where v is the velocity of the propagation, whose possible values will be 
determined later. By substituting (J3J in (J2j) and assuming s / 0, we get an 
ordinary differential equation (ODE) to be satisfied by </>(£): 

d 3 cp fid 2 (j) fv a (3 2 



d£ 3 s d£ 2 \s s s J d^ 

A further simplification of this equation is accomplished by means of the 
following linear transformation of the dependent and independent variables 

e= £ 0, m =—vV), (5) 

fi as 

which transform Eq. Q in 

d 3 w d 2 w , 2 \ dw 

^T-^T + (p-2^-3^)^ = 0, (6) 



with 

vs 4/3/i 2 

P= ^ q= 3^- (7) 

In the sequel, we will use the simplified form ©, which corresponds to the 
compound KdVB equation (J2J) if q ^ and to the usual KdVB equation (jl]) 
if q = 0. It is quite obvious that Eq. © admits a first integral 

d 2 u> dw , 9 q N 

^-- + (^-^-g W 3 )=fe, (8) 

where A; is an arbitrary constant. In the following section we will summarize 
the factorization method to find solutions of general second order differential 
equations of the type (jHJ)- Then, in sections 3 and 4, we will apply this 
method to KdVB equations. 



2. FACTORIZATION OF A SECOND ORDER 
ODE 

The factorization method is a well known technique used to find solutions 
of ODE. It goes back to some papers of Schrodinger in which he solved some 
particular examples of the equation bearing his name [T5], and it was later 
developed by Infeld and Hull 12] (for more details, see also 17 and the 
citations quoted therein). 

In a more general context, the factorization of non-linear second order 
ODE has been previously studied in ^J ^3 for equations of the form 

d 2 U dU WTTN 

___ + F((/) = , (9) 

where F(U) is a polynomial function (this is precisely the form of the two 
KdVB equations (jHJ) considered in the previous section). In these works, Eq. 
|HJ) was factorized as 

[D-f 2 (U)][D-f 1 (U)]U(6)=0, (10) 

where D = d/d9. The following grouping of terms has been proposed in |19j 

d 2 U 



d6< 



^ + IH^ + /lM, = ' (11) 



and comparing Eq. @ with Eq. (|11|). we get the conditions 

h(U)h(U) = !^P-, (12) 

hm+ wm. =1 . (13) 

Then, according to ^5j, any factorization like (|1U|) of an ODE of the form 
given in Eq. Q, allows to find a compatible first order ODE 

[D -h(U)]U = $jL -h{U)U = 0, (14) 

whose solution provides a particular solution of ©• In other words, if by 
some means we are able to find a couple of functions fi(U) and fi{U) such 
that they factorize Eq. Q in the form ()1(J|) . solving the ODE Q14JI we will 
get particular solutions of © . The advantage of the factorization presented 
here is that the two unknown functions fi(U) and fy(U) can be found easily 
by factoring the polynomial expression (|12|) in terms of linear combinations 
in rational powers of U. 

Now, let us apply this technique of finding particular solutions to the 
KdVB equations considered in the previous section. 

3. FACTORIZATION OF THE KdVB EQUATION 

Let us consider first the KdVB equation in the integrated and simplified 

form given in (jSJ 

d 2 w dw , ,> , , . 

— -- + { pw - w )=k. (15) 

It can be easily proved that this equation can be factorized in the form (|1U|) 
by means of a linear combination of powers of w only if k = 0, which is 
a very restrictive condition. To circumvent this constraint, we propose a 
simple displacement on the unknown function 

w{d) = U{6)+5, (16) 

where 5 is a constant to be determined. With this change, Eq. (|15j) becomes 

«f-f + (6,-^-^-0, (17) 

where we have imposed 

k = P 5-5 2 . (18) 



Now, we can factorize (|17|) following the procedure exposed in Section 2 with 
F(U) = (p - 25)U - U 2 . Taking into account Eqs. (|T2]i - (fT3)> . the functions 
fi(U) and f2(U) must be of the form 

AU 1/2 + B, 



h{U) 



(l-B)--AU 1 ' 2 , 



(19) 
(20) 



with the following values of the two parameters A and B 

A 2 = -, B = - 
3' 5 

together with a constraint between the parameters p and 5, introduced in 



(21) 



Eqs. (Q) and (jTB]) . respectively: 



p = 25 + 



_6_ 

25 



(22) 



As we see from (|2*2*|) . the coefficient p (and therefore, according to 0, the 
velocity v of the traveling wave) depends on the parameter 5, which is arbi- 
trary since it is fixed by the integration constant k in H18|) . 

Another important remark is that if we substitute this value of p in (|17|) , 
we can realize that we get a kind of universal equation, in which all the 
parameters have disappeared. 

Particular solutions of (|17|) are obtained by solving Eq. (|14[) in this 
particular case, which turns out to be the following couple of first order 
ODE of Bernoulli type: 



whose solutions, for the negative and positive sign, are 



U-(6) 



3 

50 



U + (9) = — 
y ' 50 



1 + tanh 



1 + coth 



0o 



10 



10 



(23) 

(24) 
(25) 



In Figs.^andElwe show a plot of these universal solutions. 

Returning to the original function w{6) in ()16|) is trivial taking into 
account (|22|) . Then, using the definitions introduced in Eqs. ©, © and 
(|7|l , we get the following family of regular solutions 



u (x,t) 



v 3fi z 
a 25as 



1 + tanh 



\x{x -vt- £q) 
10s 



1 2 



(26) 
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Figure 1: Plot of the universal regular solution U (9) given in (|24|) with 

e = o. 
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Figure 2: Plot of the universal singular solution U + (6) given in (|25|) with 
#o = 0. 



and the following family of singular solutions 

-2 . (27) 



u + (x,t) = - + 



2 ( r /../_ ..^ a w n2 



a 25as 



1 + coth ' ^ ~ Vt ~ ^ o) 



10s 



Notice that after all these steps we have obtained a couple of families (one 
regular and one singular) of particular solutions of the KdVB equation Q 
depending on two free parameters, which are the arbitrary constant £o and, 
what is more important, the velocity v of the travelling wave. In other words, 
our result indicates that for every value of the velocity we have a family of 
regular particular solutions of the KdVB equation l[T|). whose particular 
features depend on the values of the parameters of the original equation. 
However, in each case the solutions are always kink-type solutions due to the 
presence of the universal solution U~(9) given in (|24j) . which is represented 
in Fig. [I] 

It is interesting to compare our results with others obtained by means of 
more elaborated or computerized methods. For instance, in pQ the so-called 
tanh-method gives a number of particular solutions called 'regular', 'singular' 
or 'special'. However, all these solutions are obtained in our result (|26|) for 
particular choices of the parameters: their 'regular' solution (20) is obtained 
when we choose fj, 2 = 100s 2 , their 'singular' solutions correspond to the 
elections fi 2 = 25s 2 or fi 2 = 100s 2 , and the 'special' solution is found when 
the velocity of the solitary wave is v = 6// 2 /(25s). Similarly, in [2171 l2~T | I22] . 
using the 'first integral method', particular solutions are found and reported 
for specific values of the velocity. From the standpoint of these methods, 
our solutions (|26j) and (|27|). as well as other particular solutions, can be 
obtained in the standard way that makes use of the elementary symmetries 
of the KdVB equations such as shifts and Galilean transformations. On the 
other hand, from the standpoint of the factorization method the equivalent 
of these symmetries reduces to attributing relative values to the parameters 
/j, and s as stated above. 

It is worth noting that if we consider the possibility of having imaginary 
constants 9q in (|24 j) -(|25 jl or £o in (|26 [) -(|27j ) . both types of solutions would 
be in fact two versions of a unique solution. It is easy to show that such 
imaginary constants are not forbidden mathematically by the initial condi- 
tions. For example, the change 9q — ► #0 + hiix in (|24"|) produces (|25|) . More 
interesting situations are obtained for other complex values, for instance if 
we choose 9q = — bin/ 2 we get a complex solution of the KdVB equation, 
whose real and imaginary parts are represented in Figs. El and |IJ respec- 
tively. The behavior of this complex solution, having the periodicity of the 



trigonometric tangent introduced through the imaginary phase, is quite ap- 
pealing and is straightforwardly related with the complex scalar solutions 
mentioned in Refs. [21HI21] . 
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Figure 3: Plot of the real part of the universal regular solution U (9) given 
in (gU) with O = -5wr/2. 



0.15 




0.15 



0.1 



0.05 



-20 -10 10 20 30 



Figure 4: Plot of the imaginary part of the universal regular solution U (9) 
given in (jUJ) with 9 = -5i7r/2. 

In fact, if we write 9q = iair (denoting now the solutions as U~{9) = 
U~(6,a) to take into account the initial condition), and examine the para- 
metric dependence of the real and imaginary part of the solution, we see 
that the real part is passing from the pure kink type (for a = 0) to the 
singular solution for a = — 5, whereas the imaginary part evolves from the 
pure bell shape to more complex shapes at increasing a, see Figs. [S] and El 
respectively. 

Our three-dimensional plots are in agreement with Theorem 1 of Liu and 
Liu [2H] (see also their Remark 3) on complex KdVB solutions that states 
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Figure 5: Plot of the real part of the universal regular solution U (9, a) 
given in (J2U). 
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Figure 6: Plot of the imaginary part of the universal regular solution 
U~(9,a) given in l|22|) . 



that the real part is kink like whereas the imaginary part is bell-shaped. 
However, we have found here that this is a consequence of the transition of 
the regular solution towards the singular one and vice versa when varying 
the amount of radians in the imaginary phase. In addition, if one looks 
carefully to the sections of constant a, one will find that there is a range of 
a in which kink-type solutions with a pocket on the left tail and a bump on 
the right tail occur. Such kinks have not been previously reported in theory 
although it seems that they have been observed experimentally |1U| . 

4. FACTORIZATION OF THE COMPOUND KdVB 
EQUATION 

We will undertake the same program for the travelling-wave solutions of 
the compound KdVB equation (jSJ with g/0. We can factorize it directly 
as in Eq. (|1U[). taking into account that now w = U and F(U) = pU — U 2 — 
qU 3 — k, by choosing the factors in the form 

fi(U)U = AU 2 + BU + C, (28) 

f2(U) = 1 - ^^ =-2AU+(l-B). (29) 

The restriction conditions (|12 [) -(|13j ) lead us to the following values of the 
parameters: 

A 2 = §, B = ±±±, (30) 

together with 

^C = k . (32) 

The last relation (|32|) means that the coefficient C can take any value, since 
k is an arbitrary integration constant. This fact means that we can get 
factorizations for any value of p, i.e., for any velocity v of the travelling 
wave. 

Particular solutions of (jSJ are obtained by the corresponding first order 
ODE (|14j) . that in our case is 

— - AU 2 -BU-C = (33) 

da 

10 



once the parameters A, B,C are replaced by their values <|3U|) (j31|) . This 
is a simple Riccati equation whose solutions are well known. One type of 
solutions is given by 



uHe) 



-1 1 
± 



where 



Sq 3^/2q~ 



A 



1 + A tanh 



A(fl - ) 
6 



18p + 



3. 



(34) 



(35) 



Going back to the original variables of Eqs. ®, (JSJ) and JJJ), the cor- 
responding solutions of the original compound KdVB equation (J2J) can be 
written as 



u (x,t) 



a 



V 



2/3 V6^i 



1 + A tanh 



/jlA(x— vt— £o 
6s 



(36) 



A 



'I8vs 9sa 2 

—7T- + TTTT^ ~ 3. 



ft' 



2/V 



(37) 



A plot of the functions u + (x, t) is given in Fig. [7] for a certain choice of the 
parameters of the original KdVB equation and different values of the wave 
velocity. 



-10 









0.5 










-0 5 






l£^=^ 


-1 


^^f 








-1.5 









10 



Figure 7: Plot of u + (x,t) = u + (^) given in 1)36|) for the following val- 
ues of the parameters: a = 3, /3 = 2, fj, = 1, s = 2, £o = an d 
v = -1.04,-1.01,-0.94,-0.74,-0.54,-0.04, being the constant solution 
the corresponding to v = —1.04. 
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A second type of rational solutions of (|3*3")) appears when A = 0, that is, 
for the particular value of v given by 

2 a 2 a 2 

6p = 1 or v = -. (38) 

y q 6s 4/3 V ' 

This singular solution is 

± kp/A A + l 

u ° {9) --AThe~-6A^' (39) 

where A = ±-^/q/2 and k$ is an integration constant. For the special case 
where the free constant k$ = 0, we get a constant solution: 

Notice that this constant is directly obtained as a limit of the solutions (|34|) 
when A — v 0. If ko / 0, the second type solution (|40j) is recovered taking 

A 3ni , N 

"» = -^-x (41) 

and then the limit A — > 0. 

The rational solution (|39j) expressed in the initial variables becomes 

",T(.'-./)- — ( 1± \/^r| r^ 2 , (42) 

where v is given in (]38[) 3 . 

The first type solutions (|34|) were obtained in jS] (with some mistakes, 
that have been corrected here) using the 'homogeneous balance method', 
and also in other papers [U 03 ESS [^ [H] using different methods, but after 
a considerable amount of computation. However, we have not found any 
reference to the rational solutions. 




3 If we introduce the parameter e = /i. / 2/ ^ we can write the rational solution in the 
simpler form 



u (x,t) = - — 



l±e + 6eko 



e± k (x — vt — £o) 



where v= [j£\ [e 2 - 1]. 
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5. FINAL REMARKS 

In this paper, the factorization method has been shown to be well adapted 
in the search of particular solutions of non-linear differential equations of 
KdVB type. This method has a number of advantages compared to others 
applied to the same problem: i) the basic concept is quite simple and fol- 
lows the same pattern already used in linear equations; ii) the computations 
needed to develop it are quite straightforward; iii) it allows to use analyt- 
ical arguments, for example when looking for the rational solutions of the 
last section. Thus, the factorization technique is an alternative method to 
other well-settled procedures and, as we have shown here, can be success- 
fully used to get exact particular KdVB solutions. It can be applied to other 
non- linear partial differential equations, looking in an efficient way for par- 
ticular solutions of physical interest. An immediate task is the Kadomtsev- 
Petviashvili-Burgers (KPB) equation that is currently considered to provide 
a more realistic framework for the study of non-linear wave phenomena in 
the low and higher altitude auroral region |25j . 
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